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00 . Abstract 

o , 

CN . In a recent paper [Nieto M M 1996 Quantum and Semiclassical 

00 ■ Optics, 8 1061; |quant-ph/9605032|| , the one dimensional squeezed 

and harmonic oscillator time-displacement operators were reordered 
in coordinate-momentum space. In this paper, we give a general 
approach for reordering multi-dimensional exponential quadratic 
d \ operator(EQO) in coordinate-momentum space. An explicit com- 

q-i putational formula is provided and applied to the single mode and 

double-mode EQO through the squeezed operator and the time dis- 
^ ■ placement operator of the harmonic oscillator. 
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1 Introduction 



The exponential quadratic operator(EQO) plays an important role in quantum mechanics 
and quantum optics. In quantum optics, such operators occur ubiquitously in topics 
related to coherent and squeezed states. Consequently, it has always been important 
to devise and explore simplifying computational procedures for reducing these operators 
into some manageable forms. In many applications, one usually expresses these operators 
into their normal ordered forms. In a recent paper [1], it has been shown that it is also 
convenient to consider the reordering of these operators in coordinate-momentum (x-p) 
phase space as: 

exp[<5] exp[aa: 2 ] exp[(3xd] exp[7<9 2 ] 

where S, a, f3 and 7 are c-number parameters. Such reorderings, together with the follow- 
ing identities [1,2]: 

exp[cd]h(x) = h(x + c) (1) 



exp[rxd]h(x) = h(xe T ) 



(2) 



exp[cd 2 ]h(x) 



[4ttc]V2 



00 


(y - xf 


/exp 

-00 


4c 



h(y)dy, 



(3) 



facilitate the computations of the wavefunction. Moreover, as pointed out in ref[l], re- 
ordering the operators in x-p space can be applied to systems [3] with time dependent 
potentials such as: 



V(x, t) = g {2 \t)x 2 + g (i) {t)x + g w {t) 



(4) 



Reordering EQO in x-p phase space is therefore an interesting problem that deserves fur- 
ther investigations. Following Wei and Norman [4], Nietofl] has reduced the one dimen- 
sional EQO reordering problem in x-p space into the solution of four coupled first-order 
differential equations with four unknowns. However, a direct calculation formula that 
relates the EQO to its reordered form is not available. Furthermore, the results have not 
been extended to the n-dimensional case. 
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In this paper, we start within the framework of x-p space and construct a very general 
approach which is suitable to reordering arbitrary mode EQOs to its reordered form in 
x-p space. In the following section, we will outline this general approach and summarize 
the essential steps. In section 3, we show that this general approach yields an explicit 
formula for the reordering of arbitrary one-dimensional EQO. The formula is then applied 
to the one dimensional squeezed operator and time-displacement operator of the harmonic 
oscillator. The results are the same as ref[l], but unlike ref[l], we need not solve a system 
of coupled differential equations. Finally, in section 4, we consider the reordering of EQO 
in two dimensions and apply the same technique to the two-dimensional squeezed operator 
and time displacement operator of the coupled harmonic oscillator. 



2 General Approach 

We denote the n-dimensional coordinate and momentum operators as: 

x = (xi,x 2 • ■ -x n ) ; d = ip = (d 1 ,d 2 , ■ ■ -d n ) 
The commutation rule for these operators is 

[xj, dj] 5ij. 

Without any loss of generality, we shall consider the following EQO, 



U = exp 



k*M D ~ F )r~ 

2 I F D 2 \ d 



(5) 



(6) 



where D 1 , D 2 , F are n x n complex matrices and D x = D 1 and D 2 = D 2 ; the tilde sign 
denotes the transpose of a matrix. It is convenient to introduce the symmetric matrix R 



as 



Di F \ l x 

and operator A = ^(x,d)R 



. By direct calculations, we first note 



F D 2 j * zy ' ' \ d 

that if L and M are n x n complex matrices and TV is a symmetric n x n complex matrix, 
then the following identities hold: 



-xNx, dM 



\ (xNxdM - dMxNx) = -x ■ NM; 



(7a) 
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-dNd, xM 
2 



= \ (dNdxM - xMdNd) = d ■ NM 

xLd, xM\ = x ■ LM 

~xLd, dM] = -d ■ LM. 

From the above identities, one arrives at 

'A,(x,d)k\ = (x,d)RE~ 1 K 



(7b) 
(7c) 
(7d) 

(8) 



where K is an arbitrary 2n x 2n complex matrix and E denotes ^ / o ' ^ as an 

n x n identity matrix. 

Using the above formulae and commutation relations, one can recursively compute the 
following relations 

(x,d)Ri:-\ 



A, (x,d) 

A, [A, (x, d)] ] = (x, d)^- 1 ■ i?^r 1 = (x, d) (r^ 1 )' 



Applying Baker- Campbell-Hausdorff (BCH) relations, one gets 

1 



U (x, d) U- 1 = (x, d) + [A, (x, d)\ + ^ [A, [A, (x, d) 
which immediately yields 

U (x, d) £T 1 = (x, d) ■ exp (R ■ S" 1 ) 
T T 

We next denote T — | ] = exp (R ■ S" 1 ) then we have 

T21 T 2 2 



+ 



(9a) 

(9b) 
(9c) 

(10) 

(11) 



U (x, d) £T 1 = (x, d) 
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(12) 



T21 T22 

where Tu,T±2,T2i and T22 are n x n matrices. The Ty, i,j = 1,2 matrices are not 
independent. To see this, we note that the exponential matrix, exp (i?S _1 ), satisfies: 



5T 1 exp (-R^ 1 ) S = exp (-E^i?) = exp ^ R^ 1 j 



(13) 



As E 



-i 



-E, the above equation fll3|) becomes 

exp f iffi" 1 j Eexp ^i?E _1 

which, in our notation, can be recast as 

TET = E. 



Expanding and equating the entries 


in 


eq([T|), 


we get 


^22-^11 




^12^21 


= i; 






^21^11 


= T n T 2 i 






T22T12 


= T12T22 



(14) 



(15) 



(16a) 
(16b) 
(16c) 



One can then easily manipulate eq(|l^) to get the relation 

11 — 1 22 + 1 12^22 J 21 

Furthermore, by these identities, one can always have the following decomposition 



(17) 



^ii ^12 
T21 T 2 2 



I W 
/ 



e Y 0^ 
e~ Y 



I 
Z I 



with 



W = TioTo 2 -\ Z = To^TV Y = - In T 22 . 



(19) 



Let 



and 



U\ = exp 
U 2 = exp 

U 3 = exp 



2 M U 



2 My / \ d 



w° ? 

2 l ' M Z /Id 
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Using eq(ll), one has 

Ui(x,d)U^ = (x,d) 
U 2 (x,d)U^ = {x,d) 
U 3 (x,d)U^ = (x,d) 
Thus the reordered EQO, U = UiUiU^, satisfies the following relation 

if ( ,,8) &-> ={xMi w \{ e l %)(t°)= <*,a> ( : u : 12 1 . <2o) 



/ w 
o / 

e Y 0^ 
e" y 

/ 
Z I 



As shown in the Appendix, operator U~ l U' commutes with all Xi and pi so that U differs 
from U' by a c-number factor. This factor can be shown to be unity by evaluating the 
matrix element between any two states to U and U' respectively (see Appendix or ref[5] 
for details). Finally, we arrive at the following formula for reordering the EQOs in an 
n-dimensional x-p space: 



exp 



e \(trY) e -\xWx e xYd e \dZd ^1) 



In principle, one can reorder any n-dimensional EQO in x-p space through eq(pT|). 

In summary, one can compute the EQO reordering in x-p phase space according to 
the following fixed procedure: 

1. Given any EQO, one can rewrite it in the form of eq(6) to obtain the matrix R, and 
hence the matrices Di, D 2 and F . 

F —Di 

2. One then computes the exponential matrix exp(i? ■ = exp I ~ I and 



D 2 -F 



obtains the matrix 



Tii T\2 
Tii T 2 2 
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3. By eq(19), one can construct W, Z and Y explicitly. 

4. Finally using eq(21), one arrives at the reordered form. 



3 One-Dimensional Application 



We now apply the above results to one- dimensional problems and the general procedure 
simplifies considerably in this case. For one dimensional problems, we have 



U = exp 



1 



x, d) 



a c 
c b 



x 
d 



(22) 



where a, b and c are all arbitrary c-numbers. Straightforwardly, we easily get 



exp(_R • £ ) = exp 



-a 
-c 



cosh 9 + c • sinh 
b ■ sinh 

T\i T12 
T21 T22 



—a ■ sinh 
cosh 9 — c ■ sinh 



(23) 
(24) 



where 9 = \J c 2 — ab. Using eq(19), one obtains 



W 



—a 



sinh 9 ■ (cosh 9 



_1 ; Y = — lnfcosh 9 — — sinh 9}; Z = — sinh 9 ■ (cosh# 

9 9 



(25) 



Substituting eq(g^) into eq (|2T|) gives 



exp 



a c 
c b 



x 

d 



•y/ COSh t 



(26) 



exp[|| sinh ^ (cosh ^ 



f)-V]exp[- ln(cosh^ - f)x9] exp [if sinh (cosh # - fj^d 2 ] 



Eq (|26|) is an explicit formula for reordering any arbitrary one- dimensional EQO. 

To illustrate the use of eq([26|), we consider two specific examplesfl]: the time displac- 
ment operator of the harmonic oscillator and the squeezed operator in one dimension. For 
time displacement operator of the harmonic oscillator, 

T = exph|V-9 2 )] 
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Comparing this expression with eq(EBp, we get 



a = —it; b = it; c = 0; 6 = yO 2 — (—it) ■ it = it 
Using eq(|26|) , it follows 



T 



—j= exp[— - tantx 2 ] exp[— In costxd] exp[- tant<9 2 ] 
ycost 2 2 



which is just eq(44) of ref.[l]. 

The one dimensional squeezed operator is (eq(9) of ref. [1]): 

S(z) = exp[- Zl (xd + 1/2) + iz 2 (x 2 + d)/2] 

which can be rewritten as 



U = exp 



hz,d)[ %Z2 ~ Zl ) ( x 

2 \ -z x iz 2 \ d 



Comparing it with eq(12), 



a = b = iz 2 , c = —zi, 6 = J zi 2 + z 2 2 = r 



Using eq(|26|), one easily sees that 

1 



U 



cosh r + z± sinh r 



■ exp 



— - sinh r (cosh r H — ^-sinhr) 1 x 2 

2r v r ' 



(27) 



(28) 



exp 



Zl 



— ln(cosh r H sinh r)xd 



exp 



— sinh r(coshr + — sinh r) 1 d 2 

2r r 



which is just the eqs.(37) of ref.[l]. 



(29) 



4 Two- Dimensional Application 

Finally, we consider the two dimensional problem and reorder some two dimensional EQOs 
in x-p space. The two-mode squeezed operator is given by [8] 



S = exp gctia 2 — g*ai + a 2 
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Using (ai + , Oj) = -^(^i, di) ^ ^ ^ | , we can rewrite this squeezed operator 5 as 



5 = exp 



(30) 



where a 



' ^ = 73 ( j ! ™id /' is the 2 x 2 identity matrix. Here, both 



1 



i2 = tf| ° I N~ ] 



x and (9 are two-dimensional (two modes) vectors. Let R be the matrix given by 

■g*a 
ga 

g*a " 



It is easy to see that 



exp (iffi -1 ) = iV • exp ^ 



ga 



•iV" 1 



cosh \g\ ■ I + sinh |g| • cr 



2|fl| 



sinh |g| • cr 



%^ sinh |#| - a 
cosh I <7 1 • I — sinh \q\ • a 

\i>\ 2|p| i»i 



(31) 
(32) 



Following our general procedure and denoting s± as — — — sinh \g\, we get via eq(21) 

\y - I s + cosh |^| 

cosh "lfl- s + \ cosh \g\ s + 

s + cosh | g | 
cosh "l»l- s + \ cosh 



Y = - In 



cosh |g| 



101 



(33) 



-s+ cosh \g\ 



With these quantities solved, one gets from eq(ll) the x-p reordered form for the two 
modes squeezed state operator. 

For the time displacement operator of a two-dimensional coupled harmonic oscillator 
with the Hamiltonian 



H = -(xx + dd) + Axixi, — 1 < A < 1, 



we have the time displacement operator 



U = exp 



(34) 



where R 



-itM 
itl 



and M 



1 A 
A 1 



With these notations, we see that 



exp 



/ cos (ty/M) iy/M 

ism(ty/M) , , N 

v _ 1 cos (ty/M) 

V y/M V > 




(35) 



cos u sin a; 
sin a; cos a; 



where y/M — \ I and c<j = \ sin 1 A. Again using eq(9), we have 



W = iy/M tan {ty/M 



Z = -^Un(ty/M 



(36) 



Y = - In 



cos 



(*-n/m) 



From eq(21), the x-p reordered form for the two-dimensional time displacement operator 
of the coupled harmonic oscillator can thus be written down. 



Appendix 

In this appendix, we shall show that the operator U~ 1 U' commutes with the position and 
momentum operators and consequently U differs from U' by a c-number which can be 
shown to be unity. We first note that U and U' satisfy the relation 

U(x,d)U~ x = U'ix^d)^- 1 = (x,d)T. 

For the position operator, since 

U xU' 1 = U' x U'~\ 
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we have 

xU- l U' = U^U'x. 

This means UU' -1 commutes with all position operators. Similarly, one can show that 
U~ X U' commutes all momentum operators. Clearly, by Schur's lemma, one concludes that 
U^U' is proportional to unity and thus U' = c ■ U. 

We next determine the value of c. Let \f > and \g > be the eigenstate of operator 
x and d with zero eigenvalue. Clearly, < f\x = and d\g >= 0. Further, using the 
defmitition of U', one can immediately see that 

<f\U'\g> = <f\U 1 U 2 U 3 \g> 

= < f\U 2 \g > since < f\Ui — and U 3 \g >— 
= < f\exp(~d¥x)\g> 

= < f\exp(-^{xYd + tiY})\g> 

= exp(-^trlnT 22 ) < f\g> . (37) 
We denote e tRj: 1 by the following form: 

' T u {t) T 12 {t) 

T21 (t) T 2 2(t) 

and proceed to calculate the matrix element value of e tA between < /| and \g > as 

v(t) = (f\e tA \g). 

However, we note that d\g >= and < f\x = 0, so that the derivative v'(t) is given by 

v\t)= l -{f\(dD 2 d + tiF)e tA \g). (38 



In eq(|38|) , we have used the identity dFx = xFd + trF. From the transformation property 
of operator e tA in eq(12) we get the following matrix identity : 

= (f\e tA dd\g) (39) 
= (/I^^S + f^^^^t + ^^le^b) (40) 
= f 22 {t)T 12 {t)v{t) + f 22 {t){f\dde%)T 22 {t) (41) 
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Without loss of generality, one can assume that det(T 22 (t)) 7^ 0, so that 

(0f\dde tA \g) = -T 12 (t)T 22 (t)-^(t) (42) 



which leads to 



(0\dD 2 de tA \0) = -v^tilD^T^t)- 1 }. (43) 



Substituting eq(d^) into eq (|38|) we get 



v'(t) = v(t) l -ti[F - D 2 T l2 T 22 {t)-\ (44) 



On the other hand, one sees that the derivative of the T matrix is given by 

d_ ( T n (t) T 12 (t) \ = (f ~Di \ f T n (t) T 12 (t) 
dt { T 21 (t) T 22 (t) I [ D 2 -F J \ T 21 (t) T 22 (t) 



(45) 



Immediately it follows 

= D 2 T 12 (t) - FT 22 (t) (46) 
Putting eq((5|) into eq (|44]) , one sees that v(t) satisfies the differential equation 

v'(t) = -\v{t) tr[^T 22 (tr] 

which can be integrated using the condition v (0) =< f\g > to give 

v(t) = exp[-^trlnT 22 (t)] < f\g > (47) 



Comparing eq (]37|) and eq([|7]) and remembering that U — e , the value of c-number factor 



is unity so that U = U'. 
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